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DR. JOHN MANIOTE 




DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the fact of 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con^ 
nection therewith. 
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1620 USERS GROUP PROGRAM REVIEW AND EVALUATION 
(fill out in typewriter or pencil, do not use ink) 

Program No, Date 



Program Name:, 



1. Does the abstract adequately describe what the program is and what Yes No 

it does? 
Commen t 

2. Does the program do what the abstract says ? Yes No. 

C ommen t 

3. Is the Description clear, understandable, and adequate? Yes No. 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Yes No. 

Comment . 

Are the Sense Switch options adequately described (if applicable)? Yes No. 

Are the mnemonic labels identified or sufficiently understandable? Yes No_ 

Comment 



5. Does the source program compile satisfactorily (if applicable) ? Yes No_ 

Comment 

6. Does the object program run satisfactorily? Yes No 

J0 Comment 

7* Number of test cases run . Are any restrictions as to data, 

size, range, etc. covered adequately in description? Yes No_ 

C ommen t 

8. Does the Program Meet the minimal standards of the 1620 Users 

Group? Yes No. 

C omment 

9. Were all necessary parts of the program received? Yes No. 

Comment 



10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 

Please return to: Your Name 



Mr. Richard L. Pratt Company 

Data Corporation ________________________ 

7500 Old Xenia Pike Address 

Dayton, Ohio 45432 

User Group Code 

- J ©TH_3 REVIEW FORM IS PART OF THE 1620 USER GROUP ORGANIZATION'S PROGRAM 
REVIEW AND EVALUATION PROCEDURE. NONMEMBERS ARE CORDIALLY INVI TED 
TO PARTICIPATE IN THIS EVALUATION. 
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DECK KEY 

EVALUATION 

1. Determinant Source Deck 

OF A 

2. Determinant Object Deck 

DETERMINANT 



by 



Reginald T. Sterling 
Associate Professor of Mathematics 

AIR FORCE INSTITUTE OF TECHNOLOGY 
Wright-Patterson AFB, Ohio 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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ABSTRACT 



TITLE : Evaluation of a Determinant 
SUBJECT CLASSIFICATION : 5.0 

AUTHOR ; ORGANIZATION : Reginald T. Harling , Air Force Institute of 
Technology, Wright-Patterson , AFB , Ohio 

DIRECT INQUIRIES TO: 

NAME: Reginald T. Harling - Address: Math Department, 
AFIT, Wright-Patterson, AFB, Ohio 

PHONE: CL3-7H1, Extension 29115 

PURPOSE : This program will evaluate a determinant of order kO or 
less. 

METHOD : The method used is that of Chio, and is described in 

Kunz, "Numerical Analysis, " McGraw-Hill, 1957, Chap 10. 

RESTRICTIONS, RANGE : Round-off error decrease the accuracy of 

the answer as the order of the determinant 
increases. 

STORAGE REQUIREMENTS : 20K 

EQUIPMENT SPECIFICATIONS : IBM 1620, with 20, 000 digits of memory, 

card input-output. No additional features 
are needed. 



ADDITIONAL REMARKS : The source program was written in SPS lan- 
guage, having macro-ops for the input and 
output conversion of floating point numbers. 
The present program differs from 5.0.005 
principally in that, in the present program, 
all conversions of numbers, between external 
and internal forms, are done by the machine. 
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AFIT 20 E7AIUATI0H OF A l i K TEHMIBAgP R. T. Barling 

Eovember, 1962 Programmer 



Purpose : This program will evaluate a determinant of order kO or 
less. 

Equipment Weeded ; IBM 1620, 20,000 digits of memory, and input- 
out -put. Mo additional features are needed. 

Method : The method used is that of Chio, and is fully described In Kunz, 
numerical Analysis. " MeOrsv— 1111, 1957, Chapter 10. Briefly, the 
procedure is as follows: 

a) Each element of the top row of the determinant Is divided 
by its first element, A 11, so that the new A 11 is unity, (in ease 

A U • 0, two columns are interchanged to make A 11 non-zero.) 

b) A new array is formed by the rule 

A iJ A il * A iJ ^ij 

c) Column 1 and row 1 are removed from the array. What remains 
is a determinant of order one less that the original. 

d) The reduction in order is repeated until a determinant of 
order two can he evaluated by cross multiplication. 

e) If, at any stage, the top row of the determinant should con- 
tain only zero elements, the value of the determinant is at once reported 
as being zero. This situation can be recognized by the fact that the word 
ZERO, rather than a numeric representation, is typed. 

Program Language : The program was written in SPS language, having macre- 

ops for the input and output conversion of floating point numbers. 

Data Input : The first data card must carry B, the order of the determinant, 
as a two-digit fixed point number (e.g., 03, 09, 25, etc) in card columns 1 
and 2. 

The second, third, etc., data cards must carry the elements of the 
determinant, commencing with those of row 1, and continuing to lower rows. 
These numbers must be in any external floating point form acceptable to 
the UC routine. They may be presented one per card, or any number per 
card, the numbers being separated by one or more blanks or commas. It Is 
not necessary to follow the last number on a card by a record mark. If 
more than I 2 elements are presented, only the first r of them will be read. 

Operation : 

1. Set the overflow switch to PBOGRAM, others to STOP. 

2. Put the SPS object deck of this program Into the read hopper and press 
LOAD. Press READER START to read in the last two cards. 
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AT IT 20 contd. 



3. Place the data cards in the read hopper and press START and READER 
START. It will again be necessary to press READER START to read in the 
last two cards. 

k. The machine will evaluate the determinant and type out 
"DETERMIHAlfP IS 

followed by its ralue. 

5. The program may then be restarted by placing new data cards in the 
read hopper and returning to 3 above. 

6. To recommence the program before it has printed an answer press in torn, 

HIST. STOP, RESET, IISERT, RELEASE, START, 

and return to 3 above. 

Comments ; A tenth order determinant is evaluated in about 20 seconds. 
A thirtieth order determinant requires 12 minutes. Round-off error 
decreases the accuracy as the order increases. A determinant of order 
30 was computed with an accuracy of five significant figures. 
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EVALUATION 0F A DETERMINANT 



START 


TDM 


itOI ,1 


nn/itf? ic. onum nnnni 




yp 


FACT 0NE 


00414 96 0997Q 09969 




TF 




OOii.96 96 09977 09967 




RNCD 


N-1 


00438 36 02280 00500 




SF 


N-1 


00450 32 02280 00000 




M 


N,N 
Q7 


00462 23 02281 02281 




SF 


nnli7ii 39 nnnQ7 nnnnn 

UU*t/*t p*v UUUp/ uuuuu 




TC 

1 r 




nnliftA 96 1 ftJiL.ii nnnoQ 
uuHoo £o 1 CftH uuupp 




TFM 


TNT n Q 


nnixQft 16 1 Riiii7 nnfinn 

UUHpO IO IO*tH/ UUUUU 




TFM 


ARC+93 ARRAY 


one: 10 16 nnAlii o*9Lci 

UUp IV IO UUO*r 1 UtHp 1 


CAB 


RArn 


ARFA 


nnc99 a7 n99Ri nnenn 

UUDZi p/ U^tOp UUpUU 




TFM 


Pn+11 ARFA+1C7 


UUJ JH ID UUpP/ UiftU 




BD 


PH9 


nnciiA iia nncft9 nnnnn 

UUpHD Hp UvpOZ UUUUU 




on 


PQ+11,2,10 


nnccR 19 nncc7 nnnfi9 

uuppo i£ UUpP / UUUUZ. 




R 
D 


pn 


nnc7n Lq nncL.6 nnnnn 

Uvp / U Hj uu? to UUUUU 




TC 

1 r 




nncft 9 96 nn£i9 nncc7 

UUJOi to UUCI L uupp/ 




AM 


r Kl 1 to , p , I u 


nncQk 1 1 nn£i 9 nnnni 

UUppH 1 1 UUOII uuuup 


PD1 


TFM 


nn£n6 1 £ nnnnn n*nnnn 

UUOUO 1 O UUUUU UUUUU 




DC 




00617 00009 

UUOI / UUUU4 


ABC 


1 NC 


,AREA 


0061R 16 1R6C7 fioAiil 

00630 4q ir6o6 onooo 
00637 nnont; 0*0000 
00642 00005 52283 




AM 


ABC+23,10,10 


00648 11 00641 000T0 




AM 


CNT,1,10 


00660 11 18447 OOOGl 




e 


NSQ.CNT 


00672 24 18444 18447 




BE 


AAB 


00684 46 00720 01200 




BNR 


ABC, AREA 


00696 45 00618 02283 




B 


CAB 


00708 49 00522 ooooo 


AAB 


MM 


N, 10,10 

97 


00720 13 02281 000T0 




SF 


00732 32 00097 ooooo 




TF 


TN.99 


00744 26 18492 00099 




TF 


CNT.N 


00756 26 18447 02281 




SM 


CNT,2,10 


00768 12 18447 00002 




BP 


TUFF 


00780 46 01030 01100 




BZ 


EASY 


00792 46 00898 01200 


PRINT 


RCTY 




00804 34 00000 00102 




WATY 


MESS1 


00816 39 18449 00100 




0UTC 


ARE A, ARRAY, 01 608 


00828 16 18657 50851 
00840 49 18626 00000 
00847 00005 52283 
00852 00005 52451 
00857 00005 51608 
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WATY AREA 
H 

B START+12 
EASY FM ARRAY, ARRAY+30 



FM ARRAY+10 ,ARRAY+20 



FS ARRAY, ARRAY+10 



FM ARRAY, FACT 



TUFF 



PRIN 



NZ 
A22 



B24 



B25 



B 

TFM 
TFM 
BD 
C 

BE 
AM 
AM 
B 

RCTY 

WATY 

WATY 

H 

B 

SM 

BZ 

TFM 

TFM 

TFM 

TFM 

MM 

A 

A 

TFM 
TFLS 



PRINT 

BBB+1 1 , ARRAY -9 

CNT,1,10 

NZ, 

CNT ,N 

PRIN 

CNT,1 ,10 

BBB+1 1 ,10,10 

BBB 

MESS1 
MESS2 

START+12 
CNT, 1,10 
A23 

B 24+ 28, ARRAY 
B 25+ 23, ARRAY 
B25+28, ARRAY 
B26+ 23, ARRAY 
CNT, 10, 10 
B25+28,99 
B26+23,99 
CNT,1 ,10 
TEMP, 



TFLS , 



00862 

00874 

00886 

00898 

00910 

00917 

00922 

00928 

00940 

00947 

00952 

00958 

00970 

00977 

00982 

00988 

01000 

01007 

01012 

01018 

01030 

01042 

01054 

01066 

01078 

01090 

01102 

01114 

01126 

01138 

01150 

01162 

01174 

01186 

01198 

01210 

01222 

01234 

01246 

01258 

01270 

01282 

01294 

01306 

01318 

01325 

01330 

01336 



39 02283 

48 00000 

49 00414 
16 18657 
49 18546 
00005 
00005 

16 18657 
49 18546 
00005 
00005 
16 18657 
49 18526 
00005 
00005 
16 18657 
49 18546 
00005 
00005 
49 00804 
16 01065 
16 18447 
43 01186 
24 18447 
46 01126 
11 18447 

11 01065 
49 01054 
34 00000 
39 18449 
39 18479 

48 00000 

49 00414 

12 18447 
46 01522 
16 01334 
16 01359 
16 01364 
16 01389 

13 18447 
21 01364 
21 01389 
16 18447 
16 18657 
49 18586 
00005 
00005 

16 18657 



00100 
00000 
00000 
50921 
00000 
52451 
6"2481 
50951 
00000 
52461 
52471 
50981 
00000 
52451 
52461 
51011 
00000 
52451 
52279 
00000 
52442 

ooo5 1 

00000 
02281 
01200 

ooo5i 
oooTo 

00000 
00102 
00100 
00100 
00000 
00000 

00061 

01200 
52451 
52451 
52451 
52451 

000T0 

00099 
00099 
00001 
51329 
00000 

T8502 
5oooo 
51359 



B26 



B27 



A23 



A3 



B28 



B29 







01348 


49 18586 00000 






01355 


00005 5oooo 






01360 


00005 50000 


TFLS 


,TEMP 


01366 
01378 

01385 
01390 


16 18657 51389 
49 18586 00000 
00005 5oooo 
00005 T8502 


C 


CNT,N 
B27 


01396 


24 18447 02281 


BE 


01408 


46 01492 01200 


AM 


CNT, 1 ,10 
B24+28,TN 


01420 


11 18447 00091 


A 


01432 


21 01334 18492 


A 


B25+23,TN 


01444 


21 01359 18492 


A 


B25+28.TN 


01456 


21 01364 18492 


A 


B26+23,TN 


01468 


21 01389 18492 


B 


B24 


01480 


49 01306 00000 


FM 


FACT.M0NE 


01492 
01504 
01511 
01516 


16 18657 51515 
49 18546 00000 

• ✓ IV^TW WWW 

00005 52279 
00005 T85 1 2 


FM 


FACT, ARRAY 


01522 
01534 
01541 
01546 


16 18657 51545 
49 18546 00000 
00005 52279 
00005 52451 


TF 


C0N, ARRAY 


01552 


26 18522 02451 


TF 


C0N-2,ARRAY-2 
A3+23, ARRAY 


01564 


26 18520 02449 


TFM 


01576 


16 01623 92451 


TFM 


CNT, 0,10 
,C0N 


01588 


16 18447 00000 


FD 


01600 

01612 
01619 
01624 


16 1865 7 91623 
49 18566 00000 
00005 90000 
00005 T8522 


AM 


CNT, 1,10 
A3+23,10,10 


01 630 


11 18447 00091 


AM 


01642 


11 01623 000T0 


C 


CNT,N 


01 654 


24 18447 02281 


BL 


A3 


01666 


47 01600 01300 


TFM 


1,1,10 

B 29+28, ARRAY+10 


01678 


16 18524 00091 


TFM 


01690 


16 01832 92461 


TFM 


B28+28, ARRAY 


01702 


16 01802 92451 


A 


B28+28,TN 
B30+ 23, ARRAY 


01714 


21 01802 18492 


TFM 


01726 


16 01 887 92451 


A 


B30+23,TN 


01738 


21 01887 18492 


AM 


830+23,10,10 


01750 


11 01887 OOOTO 


TFM 


CNT, 1,10 


01762 


16 18447 0OO9l 


TFLS 


AREA+10, 


01774 
01786 
01793 
01798 


16 18657 61797 
49 18586 00000 
00005 62293 
00005 60000 


TFLS C0N, 


01804 


16 18657 51827 






01816 


49 18586 00000 
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FM C0N,AREA+1O 



B30 FS ,C0N 



AM B29+28.10.10 

AM B30+ 23, 10,10 

AM CNT,1,10 

C CNT,N 

BL 828 

A B28+28,TN 

TFM B29+28.ARRAY+10 

AM B30+23,10,10 

AM 1,1,10 

C I ,N 

BL B28-12 

TFM 1,1,10 

TFM B31+23, ARRAY 

TFM B3 1+28, ARRAY 

A B31+28.TN 

AM 831+28,10,10 

TFM CNT,1,10 

B31 TFLS , 





AM 


B31+23,10,10 




AM 


B31+28,10,10 




AM 


CNT.1,10 




C 


CNT.N 
B31 




BL 




AM 


B31+28,10,10 




AM 


1.1.10 




C 


I.N 

B31-12 




BL 




SM 


N,1,10 




B 


AAB 




DC 


8,10000000 


0NE 


DC 


2,01 
10 


FACT 


DS 


N 


DS 


2 


AREA 


DAS 


80 


ARRAY 


DSB 


10,1600 


NSQ 


DS 


3 
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01823 00005 T8522 
01828 0000^ 50000 
0183*+ 16 18657 51857 
01846 49 18546 00000 
01853 00005 T8522 
01858 00005 52293 
01864 16 18657 51887 
01876 49 18526 00000 
01883 00005 00000 
01888 00005 T8522 
01894 11 01832 000T0 
01906 11 01887 OOOTO 
01918 11 18447 00051 
01930 24 18447 02281 
01942 47 01774 01300 
01954 21 01802 18492 
01966 16 01832 52461 
01978 11 01887 000 To 
01990 11 18524 00051 
02002 24 18524 02281 
02014 47 01762 01300 
02026 16 18524 Q005l 
02038 16 02121 02451 
02050 16 02126 52451 
02062 21 02126 18492 
02074 11 02126 00010 
02086 16 18447 00051 
02098 16 18657 52121 
02110 49 18586 00000 
02117 00005 50000 
02122 00005 50000 
02128 11 02121 OOOTO 
02140 11 02126 OOOTO 
02152 11 18447 00051 
02164 24 18447 02281 
02176 47 02098 01300 
02188 11 02126 OOOTO 
02200 11 18524 00051 
02212 24 18524 02281 
02224 47 02086 01300 
02236 12 02281 00051 
02248 49 00720 00000 
02267 00008 
02269 00002 
02279 00010 
02281 00002 
02283 00080 
02451 00010 01600 
18444 00003 







CNT DS 3 

MESS1 DAC 15, DETERMINANT I: 
MESS2 DAC 6, ZER0@ 



TN 


DS 


3 


TEMP 


DS 


10 




DC 


8,-10000000 


M0NE 


DC 


2,01 


C0N 


DS 


10 


1 


DS 


2 




DEND 


START 



18447 
18449 
18479 
18492 
18502 
18510 
18512 
18522 
18524 
00402 



00003 
00015 
00006 
00003 
00010 
00008 
00002 
00010 
00002 



L0AD SUBR0UT I NES 
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